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Abstract 

The quantum critical Antiferromagnetic (AFM) fluctuation spectra measured by inelastic neu- 
tron scattering recently in two heavy fermion superconductors are used together with their other 
measured properties to calculate their D-wave superconducting transition temperatures T^- To 
this end, the linearized Eliashberg equations for D-wave superconductivity induced by AFM fluc- 
tuations are solved in models of fermions with various levels of nesting. The results for the ratio 
of Tc to the characteristic spin-fluctuation energy are well parametrized by a dimensionless cou- 
pling constant and the AFM correlation length. Comparing the results with experiments suggests 
that one may reasonably conclude that superconductivity in these compounds is indeed caused by 
AFM fluctuations. This conclusion is strengthened by a calculation with the same parameters of 
the measured coefficient of the normal state quantum-critical resistivity oc T^/^ characteristic of 
gaussian AFM quantum-critical fluctuations. The calculations give details of the superconducting 
coupling as a function of the correlation length and the integrated fluctuation spectra useful in 
other compounds. 

PACS numbers: 
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I. INTRODUCTION 



Many years ago, superconductivity was discovered in heavy- fermion compounds iHsj]. It 
was suggested {4] that the superconductivity was due to collective electronic fluctuations and 
not due to electron-phonon interactions. Transport properties in the superconducting state 
we.e ^alyzed to =ho„ that =upetco„ductMty was the D-wave sy„,„,etty. It wa. 

also suggested that the D-wave symmetry is promoted by Antiferromagnetic fluctuations [7|] 
with long enough correlation lengths. This promotes scattering of fermions near the fermi- 
surface predominantly through angles around ±7r/2, which is essential for superconducting 
instability in the "D-wave" channel for a suitable fermi-surface js]. The idea of long enough 
AFM correlation lengths as essential for this mechanism is supported by the fact that in 
heavy-fermions, superconductivity occurs generally in the regime near the AFM quantum 
critical point where the correlation lengths are long but the competing AFM phase has lower 
condensation energy. 

At the same time, Random phase approximation on the Hubbard model was used to 
suggest that D-wave superconductivity is promoted in such a model {9]. The properties of the 
Hubbard model have proven controversial; there are different calculations which suggest that 
the ratio of the transition temperature Tc to the typical electronic kinetic energy parameters 
t is more than O(10~^) 10| to less than 0(10^^ ) The Hubbard model was proposed 

as a sufficient model for the cuprate compounds [l2|. But this has been questioned and the 

n n 

discovery of the predicted time-reversal breaking order parameter in the underdoped 
cuprate compounds ending at a quantum critical point, on the basis of a multi-orbital model, 
and its absence in the Hubbard Model raises doubts on its validity for the cuprates. For 
pnictides, generalization of the Hubbard model to multi-orbital situations and inclusions of 
Hund's rule couplings appears essential. 

In recent years inelastic neutron scattering in the heavy fermion compounds CeCu2Si2 
17H19| and Celrlus 20[ have provided details of the AFM fluctuation spectrum in the nor- 
mal state. The primary aim of this paper is to estimate the superconducting transition 
temperature using the parameters provided by the experiments to test the idea in detail 
for a few such compounds. To do so, we solve the Eliashberg equations for d-wave super- 
conductivity using a phenomenological AFM spectral function with which the experimental 
data is in good accord. The use of the Eliashberg equations for quantitative calculations 
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in such a situation is open to question because the Migdal expansion parameter, which is 
of O(10~^) for the electron-phonon problem is of 0(1) for such compounds if one assumes 
that the scale of the AFM fluctuations extends to the order of the electronic bandwidth. 
However, when the AFM correlation length ^ is large compared to the lattice constant a or 
{2kF)~^, the scale of the AFM fluctuations is reduced correspondingly to 0{{a/^Y)t. But 
in the limit of large correlation lengths, new questions arise jsj which are not important 
in the electron-phonon problem. The most prominent among them are the role of inelastic 
scattering in depressing Tc on the one hand [2l|, and the fact that the BCS type coupling 
constant A appears to diverge when the characteristic fluctuation frequency — ?■ and the 
BCS prefactor appears to go to 0. We have tried but have been unsuccessful in sorting out 
the results for T^. in this limit analytically. An answer to these questions and the various 
considerations for Tc from AFM interactions is however possible from the numerical solution 
of the Eliashberg equations. 

Using the measured spectra of AFM fluctuations and other measured properties, we find 
from the solution of the Eliashberg equations and comparison of the calculated Tc with 
experiments that it is reasonable to conclude that AFM fluctuations are responsible for D- 
wave superconductivity in the heavy fermion compounds. Very importantly, with the same 
parameters we calculate the measured coefficient of the anomalous oc T'^/^ contribution to 
the resistivity in these compounds. 

We note here that if one adopts that the dimensionless measure Tc/Ep for how high is 
the electronic fluctuation induced superconducting, the heavy fermions may be said to do 
very well indeed. For example, in many cases, including the compounds studies here, this 
ratio is O(10~^), similar to that of the cuprates. 

Following the proposals that AFM fluctuations may also promote superconductivity in 



the cuprate compounds 2^ , there have been many discussions of the mechanism and many 



calculations based on the Eliashberg equations. A partial list includes the following 23|. 



The most complete of t 
and Lonzarich (ML) [24 



lese calculations appear to us to be those carried out by Monthoux 



25| , both for 2 and 3 dimensional models. We present below calcu- 



lations for the 2 dimensional square lattice model with a phenomenological spin-fluctuation 
spectrum, whose results no different than those of ML. The only difference in the calculations 
is that we vary the parameters in the two dimensional model so that "nesting" at the AFM 
wave-vector quantitatively changes. The amount of nesting does have a significant change 
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in the results. But more important is that now that detailed AFM fluctuation spectra is 
available, we can use the experimental parameters to test the ideas quantitatively. We also 
discuss how to put limits on the parameters used based on sum-rule for the fluctuation spec- 
tra and show that they are inter- related while ML do not. Results for the range of physical 
parameters that we flnd relevant for the heavy fermions is not available in the published 
results of ML. This has bearing also on general conditions to determine the extent to which 
AFM fluctuations give signiflcant Tc for relevant parameters in other compounds. 

This paper is organized as follows: We present in Sec. (II) the models for fermi-surface 
and for the spin-fluctuations which we have investigated using the linearized Eliashberg 
equations. We discuss there the change of effective coupling constants with the AFM cor- 
relation length using sum-rules so that the results for numerical solutions of the Eliashberg 
equations presented later are presaged. We present the results of the calculations in Sec. 
(Ill) and discuss the important conclusions immediately after the description of the Models. 
We also present, in an Appendix, the explicit derivation of the coefficient of the T^^^ resistiv- 
ity from the measured form of the AFM critical fluctuation spectra. This is used in the text 
to estimate independently the value of a coefficient A, which is important for the calculation 



of Tc. We give the parameters that have been deduced by ine^ 



the heavy fermion compounds CeCu2Si2 17h19|] and Celrlns 20| and compare the measured 



astic neutron scattering for 



Tc with the calculations. We should emphasize that such a comparison is meant to be only 
illustrative of the physical principles involved; no detailed quantitative agreement is to be 
expected, especially given that the electronic structure of these compounds is far more com- 
plicated than assumed in the models studied. However, enough details can be provided so 
that one can conclude that the idea of AFM fluctuations near the quantum critical point in 
these compounds as the source of D-wave superconductivity is well supported. For example 
using measured properties, different levels of assumed nesting in the band-structure need a 
coupling constant A between 1.5 and 3 to get the measured Tc. In this range of A and for 
the measured AFM correlation length Tc is close to being linear in A. This range of values 
is compared with the value of A ~ 1.6 needed to get the measured coefficient of the T"^/^ 
resistivity, which is relatively insensitive to nesting. One can assess the results from the fact 
that in the range of A deduced, T^ is found to be approximately linear in A. 
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II. MODELS AND RESULTS FOR % 



A. Fermi-surface 



In our calculations we will consider two types of fermi-surfaces, a free electron fermi- 
surface and the others given by the tight binding spectrum in a two dimensional square 
lattice with various magnitudes of nesting. 

= — 2t(cos k^a + cos kyO) + At cos{kxa) cos{kya) (1) 

The fermi-surface with the tight binding spectrum are shown in Fig. [T]for four values of the 
next nearest hopping t' and the AFM wave-vector. For a given AFM fluctuation spectra, 
in the model depends on how well is the fermi-surface nested with respect to the fluctuation 
spectra. The nesting in the model changes as t' increases. In particular, we will show the 
detailed result for three fermi-surfaces, the free electron fermi-surface, the fermi surface 
(FSl) with tight binding spectrum with t' = O.At and the fermi surface (FS2) with tight 
binding spectrum with t' = O.lt. Of the four Fermi-surfaces shown in Fig. [H the one with 
t' = O.At has the worst nesting and the one with t' = O.lt has the best nesting. Fig. [2] shows 
the circular fermi-surface, FSl, FS2 and the corresponding AFM wave- vectors. 

We will discuss using the results of ML together with ours, that if properly normalized 
density of states and fluctuation spectra are used, two dimensional and three dimensional 
models give similar results for provided one adjust the ratio of the region of fermi-surface 
nesting to the total fermi-surface. This is in general is always lower in three than in two 
dimensions. 



B. Dynamical Spin Susceptibility for AFM Fluctuations, Correlation length, Par- 
tial Sum rules and Coupling Constants 

The dynamical spin-susceptibility for itinerant fermions with AFM correlations may 
be usefully divided into two parts: The normal contribution of non-interacting fermions 
XoiQ, i^) and the part Xafm(^5 i^) affected by AFM correlations. The two together obey the 
total magnetic moment sum-rule. The non-interacting susceptibility can play only an in- 
significant role in promoting superconductivity jsl and should be ignored. A partial sum-rule 
on Xafm(^, i^) in terms of the ordered moment in the AFM phase can be used to relate the 
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Figure 1: The fermi-surfaces given by the tight binding spectrum with four values of the next 
nearest hopping t' with filUng=1.05 and t=0.3meV. 

integrated fluctuations to the AFM correlation length. Double counting by using the sum- 
rule on the total susceptibility for the fluctuations and yet having (free) fermions interacting 
with such spin-fluctuations is incorrect as it over-counts the total degrees of freedom. Such 
considerations have usually been blithely ignored in most of the previous phenomenological 
work on this problem. 

The fermions interact with spin-fluctuations with a phenomenological Action 

q,kk',i,a,/3,7,<5 ijJ„ 

X will be chosen to have dimensions of inverse of energy (after subsuming a factor of 4yU^ 
in its definition). So (7 is a coupling function of dimension of energy, g for heavy fermions 
is the exchange energy between the conduction electrons and the /- local moments. Its 
meaning for d-band problems is more ambiguous, and may be best inferred from independent 
experiments, for example the resitivity above T^. 



A suitable phenomenological form for t 
relations, with which experimental results 



le dynamical spin-fluctuations due to AFM cor- 



17 



18| have been fitted, but assumed isotropic in 
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Figure 2: Three types of Fermi surfaces and Q-vector. The bhie Hue shows the circular fermi- 
surface, the black line shows FSl which is given by the tight binding spectrum in two dimensional 
square lattice with filling=1.05, i=0.3 meV and t' = 0.40t. The red line shows FS2 which is also 
given by the tight binding spectrum in two dimensional square lattice but with filling=1.05, t=0.3 
meV and f = O.lOt. 



spin-space here, is 



XorAFM 



AFM 



[q-Q 



(3) 

(4) 



where Fafm is the damping rate of the fluctuations, Q is the antiferromagnetic vector. 

The correlation length ^ is related to the deviation from the QCP which is controlled by 
parameters such as pressure, doping, magnetic fleld, etc. as well as by temperature. Fafm, Q 
and ^ may all be determined from experiments. The temperature dependence of C, has been 



studied by renorma! 



ization group (RG) |26|, m 

and by the self-consistent renormalization 



(SCR) methods 28|, |29|. It depends on the nature of the QCP (for example .^"^ oc T'^/^ near 
AFM QCP in the 3 dimensional case and with the dynamical critical exponent = 2.). The 
SCR theory derived that ^ has a Curie- Weiss behavior in 2d and high temperature regime 
near the magnetic QCP, C,~'^{T) ~ Xq cc T + 6 again due to dynamical critical exponent of 
2. 
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In the present work, we focused on the two dimensional system, and ^ will be assumed 
to be consistent with the form 

(e/a)-'(T) = {ClaY' + i^. (5) 

J- AFM 

We will find however that the temperature dependence of ^ is of insignificant consequence 
in determining T^. 

The correlation length ^ has two distinct important effects for determining T^. for d-wave 
superconductors. The linearized Elisahberg equations give that kernel for Cooper pair cou- 
pling in the d-wave channel in a square lattice is proportional to the projection of 
|^(k,k')|Mk-k',a;) 

to (cosA;^. — cosA;y)(cosA;^ — cos/c^) In spin-fluctuations theories, |(7(k, k')p has only a smooth 
momentum dependence. So, Cooper-pair coupling prominently depends depends on two 
quantities : the q = (k — k') dependent part of the total spectral weight and on the corre- 
lation length S,- It is obvious that a q-independent spin - fluctuation contributes zero to the 
Cooper channel coupling in the d-wave channel and for the same total q-dependent weight, 
its magnitude decreases with the correlation length. It is not possible to be quantiative 
on the latter without detailed calculations because it also depends on the nesting in the 
band-structure near the AFM Q. The conclusions arrived in the calculations below, see Fig. 
d?! [9]) and earlier (ML) is that for the typical band-structure rapdily increases rapidly 
with increasing correlation length and saturates for ^/a > 10. 

The other effect of ^ may be discussed through the the partial sum-rule on Xa¥m{.Q^^)i 
which determines the effective coupling constant. With the assumed Lorentzian form, it is 
necessary to introduce an upper cut-off uoc in the frequencies uj up to which the fluctuations 
extend. Using such a cut-off gives 



AFM 



71^ 



^ I dulmxiq, 

1 /rAFM(e/a)~' 



TT 

Tan 

2 



1 rAFM(^/a) 



log 



1 + 



r(e/a) 



2\ -2 



(6) 



Here, '^i{Sf) afm does not include the thermal spin fluctuations. Actually, spin fluctuation 
are actually quite suppressed for u ~ Fafm and we can simply use Uc ~ Fafm in calculations 
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Figure 3: The (^/a) ^-dependence of ^^(S'j^)AFM/wc(xo/2vr) for rAFM/'^c=l-0, 2.0 and 10.0. 



of Eliashberg equations. It is important to take into account that there are four equivalent 
AFM-vector for the two dimensional problem in the paramagnetic regime of the model, 
however strongly fluctuating it may be. This has been taken into account in the sum-rule 
by multiplying 4 by Imx(g', cj). For d = 3, the number of equivalent AFM vectors is larger 
and a correspondingly larger multiplicative factor should be used. 

In the regime of very long correlation lengths, (^/a)^ ^ 1, i.e. close to the quantum- 
critical point, the sum rule simply gives 

E(^')AFM^gx° + O(a/0^ (7) 

i 

{Sf)AFU may to a first approximation be estimated from the ordered moment (S) in 
nearby AFM phase but more properly from integration of the relevant momentum and 
frequency range of the measured fluctuations. The latter requires absolute measurements 
with polarized neutrons. Fig. [3]shows the (^/a) "^-dependence of ^^(S'j^)AFM/i^c(xo/27r) for 

rAFMM=l-0, 2.0 and 10.0. 

Let us now consider the sum-rule in the opposite limit, that the correlation length is 
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small compared to the lattice constant, i.e. the system is very far from the quantum critical 
point. Then 

E(^')afm ^ ^X'lm/af + 0(e/a)1 (8) 

i 

The consideration of the correlation length on automatically appear in the numerical 
results of the Eliashberg calculations. But it is useful to consider this effect qualitatively. 
As shown below and already done in a similar way by ML, the results of the Eliashberg 
calculations may be parametrized in terms of a dimensionless "bare" coupling constant A, 

A = g'NpXo (9) 

Xo may be determined in terms of (Sf) afm and therefore (approximately) to the ordered 
moment through the sum-rule. Using that the sum-rule becomes the total moment sum-rule 
in the limit of infinite correlation length and the maximum possible ordered moment, i.e. 
that of the AFM insulator (ignoring the zero-point effects), one concludes that in the limit 
of very large correlation lengths 

Aoo=^7'iVFf (^.')max — , (10) 

where / is the fraction of the maximum possible ordered moment. It will be seen from the 
calculations below that the limit for Aoo is reached for ^/a > 10, below which there is an 
exponential fall off of Tc/Fafm; other parameters remaining the same. 

For smaller correlation- lengths. Fig. (|3]) shows how the "effective" coupling constant Acfr 
for Tc decreases with decreasing correlation length. In considering the sum-rule in relation 
to measurements of x(g, w), it should be noted that non-interacting fermions also do have 
a correlation length of (2/cf)~^, which is usually of the order of the lattice constant. Such 
incoherent fluctuations do not have any significant role in promoting superconductivity jsj. 
So if the measured correlation length is only of about a lattice constant, it follows that it is 
not meaningful to speak of collective AFM fluctuations. 

As already mentioned these considerations are important in considering the validity of 
using the AFM fluctuation theory in various different contexts. The principal point is 
that the effective coupling constant Acfr and the correlation length are not independent 
parameters, as one may think from the work of ML, where A values from about 5 to about 
50 are used in the calculations with varying correlation lengths. Actually, one obtains for 
the considerations of the sum-rules above that for spin-(l/2) problems, even the coupling 
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constant Aoo is only of 0(1), because gNp is of 0(1) and so is the ratio Fafm/i^c- An 
independent estimate of A may be obtained from the normal state properties, for example 
the coefficient of the temperature dependence of the resistivity of non-fermi-liquid form in 
the quantum critical region. Again only A of 0(1) will be found consistent. 

It is also important to note that these BCS type coupling constants do not carry infor- 
mation on the retardation effects due to the frequency dependence of the interaction; these 
as well as the effects of inelastic scattering which are particularly important for anisotropic 
superconductors are properly treated through the numerical solution of the Eliashberg equa- 
tions. 

It should also be pointed out that in some heavy fermions, the quantum-critical fluctu- 
ations do not have the functional form given by the simple RG or SCR approximations as 



above, but displays "local criticality'' 



35| as suggested for the cuprates |36j. In this paper. 



we only consider fluctuations which are well specifled by the form given above. 



III. RESISTIVITY IN THE QUANTUM-CRITICAL REGION 

An expression of the resistivity p(T) in the antiferromagnetic quantum critical region 
suitable for heavy fermions may be derived with the following formula derived from the 
Boltzmann equation. 

P'\T) = / T^,dSps, (11) 



47r3 h 3 . 

where the integration is taken over the fermi-surface, and vp the fermi velocity. This assumes 
that the actual electronic structure near the chemical potential is sufficiently complicated 
that in the temperature region of interest, vertex corrections which lead to emphasis on 
large momentum scattering for resistivity are unimportant. In that case the scattering rate 
which determines the resistivity is the same as the single-particle scattering rate. This 
is true in a multi-sheeted fermi-surface and is suitable for heavy fermions. For a single 
weakly anisotropic band, as in the cuprates, this is in general incorrect. In such a situation. 



the resistivity 



dependence [37|. 



or large AFM correlation lengths is close to the Fermi-liquid temperature 



Equation ffTTj) can also be expressed as follows. 

P-' = ^(r,),s. (12) 
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where (■ ■ ■ )fs = I ' ' ' (^^fs means the average over the fermi- surface, m* the renormal- 
ized effective mass. 

Here, can be derived from the imaginary part of the self-energy. 

n 



2rg 



-ImS(A;,£ + z5)|,^o, (13) 
where the self-energy due to the antiferromagnetic quantum fluctuations is given as follows. 

S(p,2e„) = //cbT^^ G(p - - iwm)x(g,iwm) (14) 

The result for the resistivity in the limit .^/a — ?■ oo is derived in an Appendix A. by explicitly 
calculating the self-energy given by eq. f|T^ is 



4e2 {eF/kB)^TAFM/kB 
IV. SOLUTION OF THE LINEARIZED ELIASHBERG EQUATIONS 

The superconducting transition temperature is given by the linearized version of the 
Eliashberg Equations for the normal self-energy —iuJnZ{6j:, iun) and the anomalous or pairing 
self-energy iy(6'g, zu;„). 

id ( 



[l - Z(%, iUn)] iuJn = - / ,^ ttT 1 sgn(l]„)5(\( - p, iUn - i^m), (16) 

Jfs l-^^J "^p o 

H/(0,, = - p-^.r |: - p. - (17) 

Here a;„ are the Matsubara frequencies; g is a momentum-independent coupling matrix 
element, which has already been defined , Oj: is an angle parameterizing the Fermi surface, 
N{6f:) is the density of states at angle The p-integral is over the Fermi surface, = 
dSji/ dp is the unrenormalized velocity. 



A. Results for variation of with Parameters in the Models 

Our principal general results for on the basis of solution of the linearized Eliashberg 
equations in terms of A and the parameters in Xafm{.Q,^) are given in this section. The 
numerical evaluation is done by first simplifying the Eliashberg equations f|T6l) and f|T7j) . as 
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far as possible analytically. The final expressions for the numerical evaluation, both for the 
circular Fermi-surface and the tight-binding Fermi-surfaces are given in Appendix B. 

Figure. IHshow the (^* /a) "^-dependences of T^/T afm for circular fermi- surface, and Fig. 
E] gives the magnified results in the region ^*/a < 10"^ of Fig. HJ Consistent with earlier 



discussions 



24 



25] . Tc/Fafm shows a very shallow peak at around ^ ~ 5 x 10 



^c/rAFM shows a drop-off as the correlation length decreases, while it also shows moderate 
decreases as the correlation length increases. ML pointed out that this moderate decrease is 
caused by the rapid diverges of Z as the correlation length increases. We note also that the 
quantum-classical crossover correction to the correlation length proportional to the factor 
7 in Eqs. (jS]) has a negligible effect on Tc. This will not be considered in any further 
calculations. 

The principal message from Fig. @ is that the infinite correlation length result for Tc is 
well obeyed up to ^/a ~ 10 with a very sharp fall off thereafter which will be seen later to 
be exponential. For large ^/a, no BCS type approximation for Tc is valid. The limit of very 
large correlation length is equivalent to the effective frequency of fluctuations — > 0, as may be 
seen from Eq.dl]). If we use the McMillan type approximation, in which the effective retarded 
coupling Am is proportional to the inverse of the average squared frequency of fluctuations, 
((jj^)"^, we get a divergent coupling. The Eliashberg theory yields a finite limit to Tc/Fafm, 
which depends on the bare coupling constant A. One may understand this result from the 
calculations of Allen and Dynes 38[, deduced for s-wave Eliashberg equation, that in the 
limit of a diverging coupling constant Tc oc y/Xui^^), where {u) may be taken approximately 
to be the square-root of (w^). 

The dependence on this bare coupling constant A in the large correlation length limit is 
shown in Fig. (Q. For small bare coupling A, this dependence does have the BCS form. 

Next we show the (^*/a) "^-dependence of Tc/Fafm for the four fermi surface (shown in 
Fig. [1]) in Fig. ([7]). For t' = OAt which corresponds to the worst nesting condition, Tc/Fafm 
shows results similar to that of the circular fermi- surface. As the next nearest hopping 
t' decreasing, improving the nesting condition, Tc/Fafm increases, and the shallow peak 
disappears. We can see more clearly this disappearance in Fig. ([8]) which shows the {^*/cl)~'^- 
dependence of the normalized transition temperature at (C/a)"^ = 0, Tc(^*"^)/Tc(^*"^ = 0). 



Finally, we shows the difference between circular fermi-surface and those for the tight 
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Figure 4: The transition temperature normalized by Tafm vs {^* /a)~'^ for circular fermi-surface 
and A = 1. 




binding dispersions in two dimensional square lattice. Figs. (E]), ffTU]) and fITT]) show the 
A-dependence of Tc/Fafm for circular fermi-surface, one with the worst nesting condition 
given in Fig. ([1]) which is given by tight binding dispersion, eq. ([1]) with filling=1.05, t=0.3 
meV and t' = 0.40t, one with the best nesting condition given in Fig. ([1]) which is given by 
tight binding dispersion, eq. ([1]) with filling=1.05, t=0.3 meV and t' = OAOt. As in the case 
of the /a-dependence of Tc/Fafm, the result for the worst nesting condition is similar to 
that of the circular fermi-surface. When the fermi-surface does not have a good nesting as 
shown in Fig. (fTOj) . Tc/Fafm shows behavior similar to that for the circular fermi-surface 
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Figure 6: The transition temperature normalized by Fafm vs A for circular fermi-surface and 
(r/a)-i = 0. 
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Figure 7: The transition temperature measured by Tafm vs /a)~'^ for four fermi-surfaces shown 
in Fig. ([1]) which is given by tight binding dispersion, eq. ([I]). 

although Tc/Fafm becomes smaller than the case of the circular fermi-surface. On the other 
hand, in the case with a good nesting condition as shown in Fig. (fTTl) . Tc/Fafm takes much 
larger value than the other cases. Moreover, 7 causes the decrease of Tc/Fafm in the case of 
the best nesting condition, while it causes the increase of Tc/Fafm in the case of the worst 
nesting condition. 

Figs, dl]) and (fT2l) show the (.^*/a) "^-dependence of Tc/Fafm for circular fermi-surface 
and FS2. In this case, Tc/Fafm for FS2 also takes much larger value than the circular-fermi 
surface and that with the worst nesting condition although the latter does not shown in 
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Figure 8: The transition temperature normalized by Tc{^* ^ = 0) vs ^ for four fermi- 

surfaces shown in Fig. ([TJ. 




(t/a)"' 

Figure 9: The magnified figure of Fig. [7] in the range of 10"^ < !^ 1- 



this paper. Moreover, Tc/Fafm for FS2 does not decreases even as the correlation length 
increasing compared to the other cases. Such a decrease of Tc/Fafm is caused by the excess 
increase of the renormalization Z as the correlation length increases. In other words, the 
good nesting condition suppresses the excess increase of the renormalization Z. Moreover, 
as increasing 7, Tc/Fafm decreases for all values of /a)~'^. 

ML have also presented detailed results for calculations on a 3 d electronic dispersion 
with the symmetry of a cubic lattice. They remark that other parameters being the same 
2 dispersion gives higher Tc than a three dimensional dispersion. Based on our results for 
changes in T^. in the 2 d problem, we conclude that this is because of the much better nesting 
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Figure 10: The transition temperature normalized by Tafm vs A for FSl which has the worst 
nesting condition which is given by tight binding dispersion, eq. ([T]) with filhng=1.05, t=0.3 meV 
and t' = OAOt. 




Figure 11: The transition temperature normahzed by Fafm vs A for FSl which has the best 
nesting condition which is given by tight binding dispersion, eq. ([1]) with fiUing=1.05, t=0.3 meV 
and t' = O.lOt. 

that is obtainable in model 2 d systems compared to the 3d systems for a given Q which 
spans the Fermi-surface in some (usually symmetry) direction. In fact, we can place our 
2 d-results for the very weakly nested fermi-surface over the ML results for the 3d Fermi- 
surface and find for other parameters the same that the systematics of the results for Tc as 
well as its value is very similar. 
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Figure 12: The transition temperature measured by Fafm vs ^ for FS2 which is given by 

tight binding dispersion, eq. ([I]) with fining=1.05, t=0.3 meV and t' = O.lOt. 
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Figure 13: The magnified figure of Fig. [12] in the range of 10 < {C* /a) < 10 



V. COMPARISON WITH EXPERIMENTS IN HEAVY FERMIONS 

In this section, we compare the estimates of from Eliashberg calculations with 
the experimental result in CeCu2Si2 and Celrlus. For convenience, we show the mea- 



sure intensity from Ref. (1l9l). proportional to the dynamic structure factor ^(Q,^) = 
coth(c(;/2T)Imx(Q, cj) in Fig. ( fT5l) at near Q = Qafm- 

Although the magnetic fluctuation spectrum found through inelastic scattering in 
CeCu2Si2 is well represented by the form of Eq. (jl]), the electronic structure is far more 
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Figure 14: The magnified figure of Fig. [12] in the range of 10"^ < < 1. 
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Figure 15: The measured dynamic structure factor at the antiferromagnetic Bragg vector as a 
function of frequency for various temperatures in CeCu2Si2. From Ref. (jiol ). The variation with 
Q to get the correlation length is also available in jiol ) and references therein. 



complicated than assumed here or in the 3 d calculations of ML. We have seen that 
Tc, especially in the limit of large magnetic correlation lengths depends only on gross 
parameters like A and secondarily on the amount of nesting. The comparison can only be 
very limited and can only give insight into the orders of magnitudes expected and give 
insight into the physics involved. 
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a. CeCu2Si2 [la ll9| 

To fit the phenomenological susceptibility to these results, the parameters takes the follow- 
ing values: 

Lattice constants: a,b = 4.103A; c = 9. 98 A 
AFMwave vector. Qafm= (0.22, 0.22, 1.46) 

Ordered Moment in the AFM phase : {S) ~ 0.2yUs/Ce Correlation length: ,^ ^ 25 A; 
Characteristic energy scales : 
Fafm =i 1.5 [meV]; 

Density of States: From the measured uniform (q = 0) Paramagnetic susceptibility in the 
normal state, one deduces the N{Ef) ^ (1/7) [meV]~^ /unit-cell. 
AFM spin-fluctuations parameter, xo'- 

A different expression from our Eq. (|4]) for x(g,cj) is used in Ref-jl?!) to parametrize 
the experiments for the same physical quantity. The experimental results in Ref.Jirl) are 
parametrized in terms of a quantity Xo^^ quantity F and ^. C, is the same as ours. The conver- 
sion from the quantity Xo, which we willl call Xo,s to our quantity x° is obtained by equating 



the integral over all q, u of the equation S4 in Ref.fllTl) with the integral of the same physical 
quantity given in Eq. ([6]). In the limit of (^/a)^ >> 1, one gets x° ~ Xo,s{ci/0^- The 
experimental result is Xo,s = 15.64/i^/mey. This then gives x^ ~ 0.4/i^/nnzt — cell/meV . 



We also note that the quantity F defined in Ref.fjlTI) is related to our quantity Tafm by 
^AFM = F(a)~^. 

Transition Temperature : Tc ~ 0.6 [K] 



b. Celrln^ 
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Although there is no long-range magnetic order nor any other identified broken symmetry 
competing with superconductivity in Celrlns, the resistivity of this material exhibits a non- 
Fermi liquid behavior similar to that observed in CeCoIns which locates the proximity to 
the antiferromagnetic quantum critical point. Moreover, the nuclear spin relaxation rate of 
Celrlus is also similar to CeColus. These result suggest that Celrlus is also located in the 
vicinity of a 3D antiferromagnetic quantum critical point. Recently, in this material, the 
temperature dependence of the magnetic correlation length is analyzed from the observation 
of the dynamical susceptibility by NMR experiments |20|. To fit the susceptibility to this 
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experimental result, the parameters take the following values. 

Lattice constants: a,b = 4.666A; c = 7.5168 A 
AFM wave vector. Qafm= (0.5, 0.5, 0.5) 

No antiferromagnetic long-range order has been investigated by varying the composition in 
this material. The chosen Qafm and ordered moment are taken to be that of the related 



compound CeCoIns 



32| 



Ordered Moment in the AFM phase: {S) ~ 0.15/iB 
Correlation length 

({*/a) ~ 10 at T = 1[K]|20| Characteristic energy scales of AFM: Fafm — 1.5[meV]; 
Transition Temperature: T^ = 0.4 [K] 

The experimental results show that both these materials are not far from the asymptotic 
large correlation length limit and that their Tc/Fafm are both about 0.03. For a circular 
fermi-surface, we may refer then simply to Fig. ([6]) and find that A ~ 3 gives such a result. 
For the best nested Fermi-surface, however, a A ~ 1 is sufficient as shown in Fig. ( ITTl) . 

We will estimate A in two different ways. To use the neutron scattering results, we need 
to know g besides the directly measured properties listed above. The renormalizations in the 
heavy fermion problem are such that near the critical point the AFM interaction between 
magnetic moments is of the same order as the heavy fermion bandwidth. Then g is of the 
order of the effective fermi-energy, i.e gNp ~ 1. Then we may use the experimental values 
of A^F and xo deduced from experiments above in Eq.(j9]) to get A ^ 2. 

The above manner of estimation has forced us to guess the value of g. We can estimate 
the value of A much better from the non-fermi-liquid resistivity proportional to T'^^"^ 
observed in the quantum critical regime of CeCu2Si2, whose coefficient is proportional to 
A. The resistivity p in the quantum-critical region for — )■ oo is given by Eq. ffTSj) . Using 
the values of CeCu2Si2 mentioned above, the resistivity is estimated p = 3.02 x 10~^AT'^/^ 
[flm] from eq. f fTSj) . The non-fermi liquid resistivity observed in CeCu2Si2 takes the form: 
p(T)/p300K = 0.151 + 0.071T3/2 |33| ^^leie PsooK ~ 70/if2cm js^. From the comparison of 
the coefficient of T'^/^-term in the resistivity between theoretical and experimental results, 
A is estimated as A ~ 1.6. This should be considered an important evidence for the rather 
obvious necessary idea that fluctuations that determine the normal state scattering also 
determine Tc, and of the consistency of the present calculations. The extent to which the 
calculations correctly estimate may be judged from the fact that in the range of A from 
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the different estimates for it, oc A. 
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Appendix A: Derivation of Resistivity Near the Antiferromagnetic Quantum 
Critical Point 



An expression for the resistivity under assumptions suitable for heavy fermions with a 



multi-sheeted fermi-surface and/or sufficient impurity scattering 39| is given by eq. (E 
in terms of the self-energy function Eq. ( 1T4|) . Here, we derive the relation (|T5l) explicitly. 
Substituting x(g*, iuj„i) in the spectral representation into eq. f lT4|) and carrying out the 
cjm-summation, one gets 

Taking the analytic continuation of S(p, ?£:„), the imaginary is given as 
ImS(p, e^ib) = ^ y — Imx(g, x)T\b{x - £ + ^p-_g-) (^tanh + coth 

= — — Imvfff, e — I tanh f" + coth - — --^^^^ | 
q 

= -^-J2 ^'°^AFM(g-ep-.-) _ f ^^^^ ^ ^ ^^^^ '-^r^) 

(19) 

We now consider the behavior at around the antiferromagnetic quantum critical point, i.e., 
(^/a)^^ ~ 0. In a low temperature region where the non-fermi liquid behavior appears, 
e ~ gives the dominant contribution for eq. f|T9l) . Moreover, using the following relation, 

tanh ^ - coth ^ = (20) 
2 2 smhx 
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eq. (fTOjl is transformed as 



-^^xor-M E ^-^^^^ (21) 



fc„r 



Next, we consider the g-integration in eq. ( 12T|) . Because the denominator in eq. fl2T]) has 
g'^ term, g' ~ gives the dominant contribution in the q -integration. Therefore, one gets 

-^-1 a^g'^ + sinh(ll^) 

Since the integrated function in eq. (!22l) rapidly decays as g' increases, we take Qc as 00 and 
obtain following result by easy calculation. 



8v27ra3 sinh(^) 
Substituting eqs. ( 1231) and ( |T3l) into eq. (fT2l) . p is given as 

'2v^ryp'Me2 sinh^" 
where we use na^ ~ 1 and Np = m*kp/{27i^h?). 



1 1/2 



AFMe' smh^ 



Here, we estimate the average over the fermi-surface in eq. (l24l) assuming that the fermi- 
surface is the sphere one. 



'sinhfc# 47r/c|y gj^j^ ^ 



1 1/2 



t)fs = TTIf ^^^^^Fs. (25) 



The dominant contribution in eq. ( 1251) comes from "hot" line where the relation |p| = 
|p — Ql = /cf is satisfied as shown in Fig. [T6l Assuming that the dispersion near the fermi- 
surface is given by linear dispersion, we obtain ^p^q — vpk± cos(2cr — 7r/2) = kp sin 2a, where 
k± is the deviation from the "hot" line. For one "hot" spot, the integration is estimated as 

1 , . /"'^ Vlvphi. sin 2a I 

- — -77/cFsma / ak± ; — ^ — . (2d) 

47rfc| J_,^ ^ sinh ^ ^ 



Changing the integration variable as x = wp^x sin 2a/(A;B^); eq. fl26|) is transformed as 

(A;bT)3/2 



v-p kc sin 2a 

c^x^^. (27) 



Sttw f/cf cos a Jo sinh a; 
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Figure 16: The cross-sectional circular fermi-surface and the wave vector p which satisfies the 
relation \P\ = \'P — Q\ = k-p ■ 

Now, we take the upper limit of the integration as oo because we consider the low temper- 
ature region, and eq. (p7|) can be calculated as 

(fcBT)3/^ 2v/2-l ^ /3\ ^ /3\ ^ 3(feBr)3/^ 
SvrwF^FCOsa -\/2 v2/ v2y 47rfF/i:F cosa 

Since such a "hot" point makes two rings whose total length is equal to 47r in the sphere 
fermi-surface, eq. ( l25i) is given by 



sinh%# 2eFCOsa 

Substituting eq. fl29l) into eq. f lM|) . we obtain 

3a/?, 1 3 , , 

4^2 

In this calculation, the Q- vector is given by 2A;Fsina. The Q- vector of the CeCu2Si2 is 
observed as (0.215, 0.215, 0.1458) giving \Q\ = lAda/n. Therefore, a is estimated as 
a ~ 7r/4, and we obtain the result used for the estimation of A: 

3 ah 3 /oi\ 
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Appendix B: Final Expressions for Evaluation of Tc 



Circular Fermi-surface 

For a circular Fermi-surface, it is possible to do the momentum integrals in the Eliashberg 
equations (fT6|l and (fT7|l analytically so that only a diagonalization in discrete frequency space 
needs to be done numerically. The final expressions used for numerical evaluation for the 
normal and the anomalous self-energy are: 



where 



z{e^,iwn) = 1 + 



A 



sgn(fi^ 



2n 



d% cos 2% cos 2x A- VA^ - - 



2n (52 + C2)VA2-52-C2' 



Tafm 

,2 



+ (e/a)-^ + a2(|fc|^ + |pr + |Q|^)-2a^ 



cos 



/3 = 2a'\p\^\Q\^ + \k\^-2\Q\\k\cos{e^-eQ 
A ■■ 



^-1 + + + |pf + |g|2) + 2a2|p1|g| cos(^. 



AFM 

E^ + C^ = 4a^|^|2 \\p\^ + \Q\^ + 2\p\\Q\cos{9j,-9^] 



X = tan ""^ 



\Q\ sinO^ + \p\ sm6p 



\Q\ cos 00 + \p\ cos Of! 



(32) 
(33) 
(34) 



Oq), (35) 
(36) 
(37) 
(38) 

(39) 
(40) 



Tight-Binding Fermi- surfaces 



With tight binding approximation, only some simplifications in the momentum integrals 
in the Eliashberg equations (flGll and (fTTIl can be done analytically. The final expressions 
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used in this paper for numerical evaluation are 



1 + 



u)n/ {rcT) 



^ sgn(fi^) 



A 



FS 



'-'71 "^^Tnl 



+ (e/a)-2 + a2(A:-p-Q)2' 



d'^St 

K{Un,^ni) = -2A 



1^2 (ia 



I a 



(41) 
(42) 



d^Sr, 



COS 29 j: cos 20^' 



PS {2'KYvj;jps (27r)2^;p.|Z(%,«fi„)| |^„_i]^|/rAFM + (^/a)-2 + a2(^-^-Q)2' 

(43) 

For both circular and tight binding Fermi- surfaces, the best numerical strategy to evaluate 
Tc is to cast Eqs. ( IMl) and (H3l) in the form of an eigenvalue equation for the eigenvector 

W/\u^\ 



E 



I a 



0. 



(44) 



is not Hermitian because K(ujn, ^n) includes 



It should be noted that the Matrix of eq, 
the renormalization factor Z{uJn, 0^. If the angle dependence of Z can be neglected, we can 
define K in a, form which does not include Z, and we obtain the eigenvalue equation with 
a Hermitian Matrix for the eigenvector W/\uJnZ{ujn)\- In s-wave superconductor case, such 
a situation, namely angle-independent self-energy appears. However, in the d-wave case, 
Z{iUn, 0£) strongly depends on ^g. On including Z in the kernel K, the latter is no longer 
symmetric for the frequency exchange, a;„ and ^2^- 

At high temperatures the eigenvalues of eq. (144|) are close to the negative odd integers. 
As the temperature decreases, the largest eigenvalue increases and crosses zero at transition 
temperature T = Tc. 
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